
figure;                 
f = gcf;

colormap bone;
fontname(f,'Arial');

[X,Y] = meshgrid(alphavec,betavec);
contourf(X,Y,lmat,20,'LineStyle','none'); 
hold on
plot(alphax,betax,'p', ...
    'MarkerFaceColor','white', 'MarkerEdgeColor','yellow', 'MarkerSize',13);
hold off

xlabel('Agglomeration in productivity (alpha)','fontsize',14);
ylabel('Agglomeration in amenities (beta)','fontsize',14);

ax = gca;                % define ax before using it
ax.FontSize = 14;

axis square

c = colorbar;
c.Label.String   = 'Relative Employment in Location 1 to 2';
c.Label.FontSize = 14;
